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The standard inflationary version of the origin of the cosmic structure as the result of the quantum 
fluctuations during the early universe is less than fully satisfactory as has been argued in [A. Perez, 
H. Sahlmann, and D. Sudarsky, Class. Quantum Grav., 23, 2317, (2006)]. A proposal is made there 
of a way to address the shortcomings by invoking a process similar to the collapse of the quantum 
mechanical wave function of the various modes of the inflaton field. This in turn was inspired on 
the ideas of R. Penrose about the role that quantum gravity might play in bringing about such 
breakdown of the standard unitary evolution of quantum mechanics. In this paper we study in 
some detail the two schemes of collapse considered in the original work together with an alternative 
scheme, which can be considered as "more natural" than the former two. The new scheme, assumes 
that the collapse follows the correlations indicated in the Wigner functional of the initial state. 
We end with considerations regarding the degree to which the various schemes can be expected to 
produce a spectrum that resembles the observed one. 
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I. INTRODUCTION 

In recent years, there have been spectacular advances in 
physical cosmology, resulting from remarkable increase in 
the accuracy of the observational techniques and exem- 
plified by the Supernova Surveys [l[ , the studies of large 
scale structure [2| and the highly accurate observations 
from various recent studies in particular those of Wilkin- 
son Microwave Anisotropy Probe (WMAP) These ob- 
servations have strengthened the theoretical status of the 
Inflationary scenarios among cosmologists. 

We should note however that while much of the focus of 
the research in Inflation has been directed towards the 
elucidation of the exact form of the inflationary model 
(i.e. the number of fields, the form of the potential, and 
the occurrence of non-minimal couplings to gravity to 
name a few), much less attention has been given to the 
questions of principle, as how the initial conditions are 
determined, what accounts for the low entropy of the 
initial state, and how exactly does the universe transit 
from a homogeneous and isotropic stage to one where 
the quantum uncertainties become actual inhomogeneous 
fluctuations. There are of course several works in which 
this issues are addressed 0, HI but as explained in 0, @, [1] 
the fully satisfactory account of the last of them seems 
to require something beyond the current understanding 
of the laws of physics. The point is that the predictions 
of inflation in this regard can not be fully justified in 
any known and satisfactory intcrprctational scheme for 
quantum physics. The Copenhagen interpretation, for in- 
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stance, is inapplicable in that case, due to the fact that 
we, the observers, are part of the system, and to make 
things even worse we are in fact part of the outcome of the 
process we wish to understand, Galaxies, Stars, planets 
and living creatures being impossible in a homogeneous 
and isotropic universe |l6|. The arguments and counter- 
arguments that have arisen in regard to this aspect of 
the article mentioned above have been discussed in vari- 
ous other places by now and we point the reader who is 
interested in that debate to that literature [1, Q . In the 
present work we will focus on a more detailed study of 
the collapse schemes and on the traces they might leave 
on the observational data. Nevertheless, and in order to 
make the article self contained, we will briefly review the 
motivation and line of approach described in detail in [|| . 
To clarify where lies the problem, and the way in which 
it is addressed in [|| we will review in a nutshell, the 
standard explanation of the origin of the seeds of cosmic 
structure in the Inflationary paradigm: 

• One starts with an homogeneous and isotropic 
spacetime[T3] The inflaton field is the dominant 
matter in this spacetime, and it is in its vacuum 
quantum state, which is homogeneous and isotropic 
too. The field is in fact described in terms of its ex- 
pectation value represented as a scalar field which 
depends only on cosmic time but not on the spa- 
tial coordinates, 4>o and a quantum or "fluctuating 
part" 5(j) which is in the adiabatic vacuum state, 
which is an homogeneous and isotropic state (some- 
thing that can be easily verified by applying the 
generators of rotations or translations to the state). 

• The quantum "fluctuations" of the inflaton acts as 
perturbations [l8| of the inflaton field and through 
the Einstein Field Equations (EFE) as perturba- 



2 



tions of the metric. 

• As inflation continues the physical wave length 
of the various modes of the inflaton field become 
larger than the Hubble Radius (Horizon-crossing 
as referred commonly in the literature), and the 
quantum amplitudes of the modes freeze. At that 
moment one starts regarding such modes as actual 
waves of in a classical field. Later on, after inflation 
ends, and as the Hubble Radius grows, the fluctua- 
tions "re-enter the Horizon" , transforming at that 
point into the seeds of the cosmic structure. 

The last step is usually refereed as the quantum to clas- 
sical transition. There are of course several schools of 
though about the way one must consider such transi- 
tion: from those using the established physical paradigms 
[H HI j to views advocating a certain generalization of the 
standard formalisms [T^- The two works focus con- 

cretely in a full blown quantum cosmology, and its inter- 
pretational problems, which are even more severe than 
the ones we are dealing with here. In [|| it was argued 
that such schemes are insufficient, in particular if one 
expects cosmology to provide a time evolution account 
starting from the totally symmetric state to an inhomo- 
geneous and anisotropic Universe in which creatures such 
as humans might eventually arise. 

The view taken in Q (and in this work) intends to be 
faithful to the notion that physics is always quantum me- 
chanical, and that the only role for a classical description 
is that of an approximation where the uncertainties in the 
state of the system are negligibly small and one can take 
the expectation values as fair description of the aspects 
of the state one is interested on. However one must keep 
in mind that behind any classical approximation there 
should always lie a full quantum description, and thus one 
should reject any scheme in which the classical descrip- 
tion of the universe is inhomogeneous and anisotropic but 
in which the quantum mechanical description persists in 
associating to the universe an homogeneous and isotropic 
state at all times. Thus in Q, one introduces a new in- 
gredient to the inflationary account of the origin of the 
seeds of cosmic structure: the self induced collapse hy- 
pothesis. I.e. one considers a specific scheme by which a 
self induced collapse of the wave function is taken as the 
mechanism by which inhomogeneities and anisotropics 
arise in each particular scale. This work was inspired in 
early ideas by Penrose [ll[ which regard that the col- 
lapse of the wave function as an actual physical process 
(instead of just an artifact of our description of physics) 
and which is assumed to be caused somehow by quan- 
tum aspects of gravitation. We will not recapitulate the 
motivations and discussion of the original proposal and 
instead refer to the reader to the above mentioned works. 
The way we treat the transition of our system from a 
state that is homogeneous and isotropic to one that is 
not, is to assume that at certain cosmic time, something 
induces a jump in a state describing a particular mode of 
the quantum field, in a manner that would be similar to 



the standard quantum mechanical collapse of the wave 
function associated with a measurement, but with the 
difference that in our scheme no external measuring de- 
vice or observer is called upon as "triggering" that jump 
(it is worthwhile recalling that nothing of that sort exist, 
in the situation at hand, to play such role). 

The main aim of this article is to compare the results 
that emerge from the collapse schemes considered in Q 
with an alternative scheme of collapse that can be said to 
more natural than the previous two. In this new scheme 
we take into account the correlations in the quantum 
state of the system before the collapse for the values of 
field and conjugate momentum variables as indicated by 
the Wigner functional analysis of the pre-collapse state. 

This article is organized as follows: In the Section HT1 we 
review the formalism used in analysing the collapse pro- 
cess. Section [IIII review how to obtain the wave function 
for the field from its Fock space description, which is then 
used in evaluating the Wigner function for the state, and 
the state that results after the collapse. Section IIVI de- 
scribes the details of the spectrum of cosmic fluctuations, 
resulting from such collapse, and finally, in Section fVl we 
discuss these results and those of other collapse schemes 
vis a vie the empirical data. 



II. THE FORMALISM 

The starting point is the action of a scalar field with 
minimal coupling to the gravity sector: 



S[(j>,g a b 




v a 4V b H ab -v (1) 



One splits the corresponding fields into their homoge- 
neous ( "background" ) part and the perturbations ( "fluc- 
tuation"), so the metric and the scalar field are written 
as : g = go + Sg and = cf>o + 8<p. With the appropriate 
choice of gauge (conformal Newton gauge) and ignoring 
the vector and tensor part of the metric perturbations, 
the space-time metric is described by. 

ds 2 = a( ? y) 2 [-(1 + 2^)dn 2 + (1 - 2%)8 ij dx l dx j ] , (2) 

where ^ is called the Newtonian potential. One then con- 
siders the EFEs to zeroth and first order. The zeroth or- 
der gives rise the standard solutions in the inflationary 
stage, where a{rj) = — 757^, with H] ~ (8tt/3)GV with 
the scalar potential, O in slow-regime so (f)' ~ — -^jjj^', 
and the first order EFE reduce to an equation relating 
the gravitational perturbation and the perturbation of 
the field 



V 2 * = AnG(/)' S^ = si 



(3) 
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with s = 4ttG(/)' q . The next step involves quantazing the 
fluctuating part of inflaton field. In fact it is convenient 
to work with the rescaled field y = a8<f>. In order to avoid 
infrared problems we consider restriction of the system 
to a box of side L, where we impose, as usual, periodic 
boundary conditions. We thus write the fields as 

V(V, e lfe -%(77), 77(77, x) = i e tk ' S *k(v), 

k fc 

(4) 

where % k is the canonical momentum of the scaled field. 
The wave vectors satisfy kiL = 2imi with i — 1,2,3, 
and y k (rj) = y k (f])d k + yk{rj)a\, and 7^.(77) = g k {^)d k + 
n k (rf)a k . The functions y k (rj), g k (ri) reflect our election of 
the vacuum state, the so called Bunch-Davies vacuum: 

(5) 

The vacuum state is defined by the condition a k \0) = 
for all k, and is homogeneous and isotropic at all scales. 
As indicated before, according to the proposal, the self- 
induced collapse operates in close analogy with a "mea- 
surement" in the quantum-mechanical sense, and as- 
sumes that at a certain time i}% the part of the state 
that describes the mode k jumps to an new state, which 
is no longer homogeneous and isotropic. To proceed to 
the detailed description of this process, one decomposes 
the fields into their hermitian parts as follows y k = 
i) k (v) + iVkiv), and 7r fe = 7^(77) + ^£(77). 
We note that the vacuum state |0) is characterized in 
part by the following: its expectation values (yf*' 1 ^)) = 
(jt k ' (77)) = and its uncertainties are Ay k ' = 
l/2\y k \ 2 (hL 3 ) and Att£' 1 = l/2| 5fc | 2 (^ 3 ). 
For an arbitrarily given state of the field |f2), we introduce 
the quantity d k = (Q\a k J \Q) = \d k a \e iak so that, 

(y£' I ) = V2Zl(y k d«' 1 ), {^ I )=V2^{g k df I ) 1 (6) 

which shows that it specifies the main quantity of interest 
in characterizing the state of the field. 
It is convenient for future use to define the following 
phases, f3 k = arg(?/ fe ) and 7 fe = arg(g fe ), keeping in mind 
that they depend on the conformal time 77. 
The analysis now calls for the specification of the scheme 
of colla pse determining the state of the field after the 
collapse [20j, which is the main purpose of the next sec- 
tion. With such collapse scheme at hand one then pro- 
ceeds to evaluate the perturbed metric using a semi- 
classical description of gravitation in interaction with 
quantum fields as reflected in the semi-classical EFE's: 
Gab = 8irG{T ao ). To lowest order this set of equations 
reduces to 

V 2 * fe = 8{6fa)n, (7) 
where {5<j)' k )n is the expectation value of the momentum 



field Scj)' k = Tr k /a(rj) on the state characterizing the 
quantum part of the inflaton field. It is worthwhile em- 
phasizing that before the collapse has occurred there are 
not metric perturbations [2~i| . i.e. the r.h.s. of the last 
equation is zero, so, it is only after the collapse that 
the gravitational perturbations appear, i.e. the collapse 
of each mode represents the onset of the inhomogeneity 
and anisotropy at the scale represented by the mode. An- 
other point we must stress is that, after the collapse, and 
in fact at all times, our Universe would be defined by a 
single state |fJ), and not by an ensemble of states. The 
statistical aspects arise once we note that we do not mea- 
sure directly and separately each the modes with specific 
values of k, but rather the aggregate contribution of all 
such modes to the spherical harmonic decomposition of 
the temperature fluctuations of the celestial sphere (see 
below) . 

To make contact with the observations we note that the 
quantity that is experimentally measured (for instance by 
WMAP) is AT/T(8,(p), which is expressed in terms of 
its spherical harmonic decomposition ^2 lm cti rn Yi rn (6 1 ip). 
The contact with the theoretical calculations is made 
trough the theoretical estimation most likely value of the 
«/m's, which are expressed in terms of the Newtonian po- 
tential on the 2-sphere corresponding to the intersection 
of our past light cone with the of last scattering surface 
(LSS): ^(t] D ,xd), ai m = J $(r)d, x D )Y l * m d 2 VL. We must 
then consider the expression for the Newtonian Potential 
([7]) at those points: 

^' f )=Ef|M^' (8) 
fc 

where we have introduced the factor T(k) to represent 
the physics effects of the period between reheating and 
decoupling. 

Writing the coordinates of the points of inter- 
est on the surface of last scattering as x = 
i?£>(sin#sin0, sin cos 0, cos 6*), where Rd is the comov- 
ing radius of that surface and 8, 4> are the standard spher- 
ical coordinates of the sphere, and using standard results 
connecting Fourier and spherical expansions we obtain 

a ^ = E "-^fj J (Wk)e rk - s Yi m (e, 4>)d 2 n (9) 

k 

= ^E^?^) 4 ^(i^^) r ^^)- ( 10 ) 

fc 

As indicated above statistical considerations arise when 
noting that the equation ([8]) indicates that the quantity 
of interest is in fact the result of a large number (actually 
infinite) of harmonic oscillator each one contributing with 
a complex number to the sum, leading to what is in effect 
a two dimensional random walk whose total displacement 
corresponds to the observational quantity. Note that this 
part of the analysis is substantially different from the 
corresponding one in the standard approach. In order to 
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obtain a prediction, we need to find the most likely value 
of the magnitude of such total displacement. 
Thus we must concern ourselves with: 



\ai 7 



16s 2 tt 2 ^ T(fc) T(fc') 



11 L 6 ^ k 2 k' 2 

kk' 

{5^' k ){5^ kl Yj l {kR D ) 3l (k'R D )YUk)Y lm {k'), (11) 



and to obtain the "most likely" value for this quantity. 
This we do with the help of the imaginary ensemble 
of universes [111 and the identification of the most likely 
value with the ensemble mean value. 
As we will see, the ensemble mean value of the product 
{5(f>k){5(f>k')* , evaluated in the post-collapse states [23^ . 
results in a form KC(fc)<5gg, , where k = HL 3 k/(Aa 2 ) and 
C(k) is an adimensional function of k which codifies the 
traces of detailed aspects of the collapse scheme. We are 
thus lead to the following expression for the most likely 
(ML) value of the quantity of interest: 



I |2 i 
\ a lm\ M L = s 



47T 2 fi ^ C{k)T{kf . 



E 



L3 a 2 ^ fc 3 
k 



jf{\k\R D )\Y lm {k)Y 



(12) 

Writing the sum as an integral (using the fact that the 
allowed values of the components of k are separated by 
Ak, = 2n/L): 



\ a lm\ 2 ML 



h f C(k)T(k) 2 . 



2na 2 



fc 3 



jfdklR^lYUk^d'k 



(13) 

The last expression can be made more useful by changing 
the variables of integration to x = kRrj, leading to 



s 2 h fC(x/R D ), 



\ a lm\ML 



2ira 2 



T(x/R D ) 2 jf(x)dx. (14) 



With this expression at hand we can compare the expec- 
tations from each of the schemes of collapse against the 
observations. We note, in considering the last equation, 
that the standard form of the spectrum corresponds to 
replacing the function C by a constant. In fact if one 
replaces C by 1 and one further takes the function T 
which encodes the late time physics including the plasma 
oscillations which are responsible for the famous acous- 
tic peaks, and substitutes it by a constant, one obtains 
the characteristic signature of a scale invariant spectrum: 

I 1 2 1 

\ a lm.\ML 00 1(1+1) ■ 

In the remaining of the paper we will focus on the ef- 
fects that a nontrivial form of the function C has on the 
predicted form of the observational spectrum. 



III. PROPOSAL OF COLLAPSE a la WIGNER 

As indicated in the introduction, the schemes of col- 
lapse considered in the first work following the present 
approach, Q, essentially ignored the correlations be- 



tween the canonical variables that are present in the pre- 
collapse vacuum state. In the present analysis, we will 
focus on this feature, characterising such correlations via 
the Wigner distribution function [12j , and requiring the 
collapse state to reflect those aspects. The choice of the 
Wigner distribution function to describe these correla- 
tions in this setting is justified by some of its standard 
properties regarding the " classical limit" (see for instance 
j!3j). and, by the fact that there is a precise sense in 
which it is known to encodes the correlations in question 
[3] ■ The Wigner distribution function for pure quantum 
states characterized by a position space wave function 
^{q) is denned as: 

W{q,p) ^ihl rf y**(9 + y)*(?-^) ex p(x)' 

(15) 

with (q,p) corresponding to the canonical conjugate vari- 
ables. 



In our case the wave function for each mode of the field 
(characterized by its wave vector number fc) corresponds, 
initially, to the ground state of an harmonic oscillator. 
It is a well known result that the Wigner distribution 
function gives for a quantum harmonic oscillator in its 
vacuum state a bi dimensional Gaussian function. This 
fact will be used to model of the result of collapse of 
the quantum field state. The assumption will be that 
at a certain (conformal) time r] k the part of the state 
characterizing the mode fc, will collapse (in a way that 
is similar to what in the Copenhagen interpretation is 
associated with a measurement), leading to a new state 
| Q) in which the fields (expressed by its hermitian parts) 
will have expectation values given by 



] K k cosQ k , (7Tfc" r )o = xj^'^Afcfcsinefc, 

(16) 

where x^' 1 * 1 is a random variable, characterized by a 
Gaussian distribution centered at zero with a spread one; 
A k is given by the major semi-axis of the ellipse char- 
acterizing the bi dimensional Gaussian function (the el- 
lipse corresponds to the boundary of the region in "phase 
space" where the Wigner function has a magnitude larger 
than 1/2 its maximum value), and & k is the angle be- 
tween that axis and the y k yI axis. 



Comparing (O with (TIB)) we obtain, 



\d k J \cos{a k + (3 k ) 



1 R,I 

V2\Vk\ Xk 



AfeCosBfe, (17) 



\d k J \cos(a k +7 fe ) = — 2- — -x k J A k ksm<d k . (18) 



V2\ 



9k\ 



From this expressions we can solve for the constants 
dk' = \d k ,I \e lah - In fact using the polar representation 
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of the y k and g k we find 
tan(afc — krj) 

obtaining 



k 2 r] c \y k \ smQ k 



g k \ co&Q k yJl + W^fY - k\y k \ sin6 fc 

(19) 



This has the form of a bi dimensional Gaussian distri- 
bution as expected from the form of the vacuum state. 
The cross term is telling us that the support of Wigner 
function is rotated respect the original axes. Rescaling 
the 7Tfe-axe to lift = ir k /k and doing a simple 2D ro- 



iR,I 



tation (i.e. y k 
tRJ o RJ 



R,I r\ 

y k ' cosOfc 



nR,I ■ 

IV sinBfe, 



n 



iR.I 



x^A k V 1 + k Wc 



k 3111 "-k - 

cosO/j — y 1 ^' 1 sinOfc) we find the principal axes of 



krj c 



\y k \ 2 k 2 sin 2 e fc + |g fc | 2 cos 2 e fe 



2|y fe ||5fe|fccose fe sin0 fc 



the Wigner function 



,RJ u ,RJ , v ) = 2(1 



1/4 



(1 + fc 2 7 ? 2)l/2 



k 2 r\ 2 



(20) 



where in all of the expressions above the conformal time 
r] is set to the time of collapse 77* of the corresponding 
mode. 

In order to obtain the expression for A k it is necessary 
to find the wave-function representation of the vacuum 

R I 

state for the variable y k ' . Following a standard proce- 
dure, we apply the annihilation operator, a fl ' 7 , to the 
vacuum state |0), obtaining the well-known equation for 
the harmonic oscillator in the vacuum state, and from 
the result we extract the wave function of the /c— mode 
of the inflaton field: 



R,I 



R,I 

v k ,v 



2k 



\ 



1/4 



I 

krj 



irhL 3 



cxp 



\ 



HL 3 {l + — 
krj 



R,i 
Vk 



(21) 



/ 



We next substitute this in the expression for the Wigner 
function, W(y^' 7 , "ftj?' 1 , 77), obtaining, 



2 1 



k 2 r\ 2 



1/4 



exp 



(— 

\krjhL 3 



Vk 



R,I R,I 



cxp 



exp 



2fc ( R,I 
hL 3 ( V * 

(1 + k 2 r) 2 ) 
2hL 3 k 3 ri 2 



(22) 



cxp 



,R,I 

Vk 



Vk 



ir 



R.I 



exp 



am 



(23) 



with the corresponding widths given by: 
4hL 3 hf 



1 + hk 2 r\ 2 



CT n' = 



a/I + lOfcV 
4hL 3 krj 2 



9&V 



1 + 5fcV - V 1 + + 9fcV 



(24) 



(25) 



Note that <7ri' > oy . The rotation angle, 9 k is given by 



26 A 



arctan 



Akrj 



1 - 3k 2 ?j 2 



(26) 



It is clear then that A k = 2cr n < 



Substituting Tt k in 8<f>' k (defined by the equation (|8|)) and 
calculating the expectation value of it in the post-collapse 
state, |f2), we obtain 



fe/n 



1 r 
a 



,R \cos (a R 



k 1 Lun V"fc t Ik + A fc ) 
i|4|cos (a£+ 7fc + A fc ) 



(27) 



where we have defined the "collapse to observation delay" 
from the collapse time of the mode k, rj k as A*, = k(r)—i] k ) 
where 77 represents the time of interest which in our case 
will be the "observation time" . 

Inserting the equation (|20p in the last expression, we can 
rewrite (S(f>' k )n as 



6 



fe/n 



krj c \/tiL 3 k 



x k + ix l 



iVc) (i + lofc 2 ?? 2 + 9fc 4 ^)V4 r + 5fc2 ^ _ ^ /1 + 10jfc2j? 2 + gkijl i 

^k\J V 1 + lOfc 2 ?? 2 + 9k 4 T]J- 1 + 3fc 2 7y2 + 



cos 



sin At 



yV 1 + 10fc 2 ?y 2 + 9fc 4 r? 4 + i _ 3^2- 



— J v 71 + 10fc 2 ?7 2 + 9fc 4 7? 4 - 1 + 3fc 2 7y 2 
kr) c 



(28) 



Now we take the ensemble mean value of the square of 
(<5</4)n, taking out a factor of k (remember that k = 
hLPk/Aa 2 , see last section) and call it C W i gner {k) 



Cyjigner (^) 



32z 2 



v/1 + 165= + 9z 4 
1 



1 + 5z 2 - v/1 + 10z 2 + 9z 4 



1 + 10z 2 + 9z^ - 1 + 3z^ 



cos At 



sin At 



■-k 



sin z A fc 



1 + 10z 2 + 9z 4 



3z 2 - 7 



8zfe cos Afc sin Aa 



(29) 



where we replaced fcrfu(fc) by z k . Henceforth (TH)) is 



\0-lm\ML 



s 2 h 

2na 2 



a 



wigner 



[x/Ri 



-T(x/R D ) 2 jf(x)dx. 

(30) 

Now we are prepared to compare the predictions of the 
various schemes of collapse with observations. 
Before doing so it is worth recalling that the standard 
results are obtained if the function C is a constant, and 
to mention that it turns out that in order to obtain a 
constant C (in this and any collapse scheme) there seems 
to be a single simple option: That the z k be essentially 
independent of k indicating that the time of collapse for 
the mode k, rj k should depends of the mode frequency 
according to T}% = z/k. For a more detailed treatment we 
refer to the article 



IV. COMPARING WITH OBSERVATIONS 

This is going to be a rather preliminary analysis concen- 
trating on the main features of the resulting spectrum 



and ignoring the late time physics corresponding to the 
effects of reheating and acoustic oscillations (represented 
by T(k)). Actual comparison with empirical data requires 
a more involved analysis which is well outside the scope 
of the present paper. 

We remind the reader that C(k) encapsulates all the im- 
print of the details of the collapse scheme on the obser- 
vational power spectrum. 

The functional form of this quantity for the scheme con- 



sidered in this article, C W i gner (|29|). has a more com- 
plicated form than the corresponding quantities that re- 
sulted from the schemes of collapse considered in [f| . Here 
we reproduce those expressions for comparison with the 
scheme considered here and with observations. In the first 
collapse scheme (f3"Tj) , the expectation values for the field 
yk and its canonical conjugate momentum ttu after the 
collapse are randomly distributed within the respective 
ranges of uncertainties in the pre-collapscd state, and arc 
uncorrelated. The resulting power spectrum has 

Ci(fc) = l + ^sin 2 A fe + — sin(2A fe ). (31) 

z k z k 

The second scheme considered in [(| only the conjugate 
momentum changes its expectation value from zero to a 
value in such range, this second scheme is proposed since 
in the first-order equation ([7]) only this variable appears 
as a source. This leads to a spectrum with 



C 2 (k) = 1 + sin" A fe 1- 



1 



-— sin(2A fc ). (32) 

Zk 



Despite the fact that the expression for C W i gner looks by 
far more complicated that C2, their dependence in z k is 
very similar, except for the amplitude of the oscillations 
(see figures 1(b) and 1(c)). Another interesting fact that 



can be easily detected in the behaviour of the different 
schemes of collapse is that if we consider the limit z k — > 
±00, then Ci(fc) — * 1 and we recover the standard scale 
invariant spectrum. This does not happen with C2(k) or 
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(a)Ci, the two field variables (j/j.) and (7Tfc), collapses 
to a random value of the dispersion of the vaccum 
state independently 




(b)C2, this scheme is proposed taking in account the 
fact that only (7Tfc) appears in the EFE at first order. 




{c)C w i gner , this scheme proposes a kind of correlation 
between the post-collapse values taking the Wigner 
functional of the vaccum state as an indicator of this 
correlation. 

FIG. 1: Plots of the three collapse schemes, we could apreciate 
that C2 (middle) and C W i gne r have a similar behavior despite 
their dissimilar functional form. 



C wi gner(k) (SCO figure [TJ . 

We recall that the standard form of the predicted spec- 
trum is recovered by taking C(k) = 1. Therefore, we can 
consider the issue of how the various collapse schemes ap- 
proach the standard answer (given the fact that the stan- 
dard answer seems to fit the observations rather well). In 
particular we want to investigate how sensitive are the 



predictions for the various schemes, to small departures 
from the case where is independent of k, which as we 
argued above would lead to a precise agreement with the 
standard spectral form. In order to carry out this anal- 
ysis, we must obtain the integrals (fT4|) for the various 
collapse schemes characterized by the various functions 
Ci(fc),C2(fc) and C W i gner (k). It is convenient to dchne 
the adimensional quantity z x = xN(x), where x = kRo 
and N(x) = TI%r x \/RD- We will be working under the 
following assumptions: (1) The changes in scale during 
the time elapsed from the collapse to the end of inflation 
are much more significant than those associated the time 
elapsed from the end of inflation to our days, thus we will 
use the approximation = —z x ; (2) We will explore the 
sensitivity for small deviations of the "zk independent of 
k recipe" by considering a linear departure from the k 
independent Zk characterized by z x as z x = A + Bx in 
order to examine the robustness of the collapse scheme in 
predicting the standard spectrum. We note that A and 
B are adimensional. 

In the figures [2 [3] and [4] reflect the way the spectrum 
behaves as a function of I , were we must recall that 
standard prediction (ignoring the late physics input of 
plasma oscillations etc) is a horizontal line. Those graphs 
represent various values of A and B chosen to sample a 
relatively ample domain. The graphs ([5l and [7]) show 
the form of the spectrum for various choices for the value 
of B keeping the value of A fixed. 

It is important at this point remind the reader -in the or- 
der to avoid possible misinterpretations- that this graphs 
are ignoring the effect of late physics phenomena (plasma 
oscillations, etc.). Our aim, at this stage is to compare 
this graphs with the scale-invariant spectrum predicted 
by standard inflationary scenarios (i.e. a constant value 
for 21(1 + l)|cv; m | 2 ) and not -directly- with the observed 
spectrum. 

As we observed before the behavior of C2 and C, 



wigner 



IS 



qualitatively similar, the main difference comes from the 
amplitude of the oscillations of the functional. 
From these results we can obtain some reasonable con- 
strains on the values of the A and B for the different 
schemes of collapse. We start by defining for a given pre- 
dicted spectrum the degree of deviation from the flat 
spectrum to be simply A imax = (j^Ej=i mox [(i(Z + 
1)^1 £m|«Z,m| 2 -S] 2 ) 1 / 2 /S where S represents the flat 
spectrum that would best approximate the corresponding 
imaginary data and is given by S = j^^'E l l z z l [ nax (1(1 + 
1) 2fp['E>m\&i,m\ 2 )- If we set a bound on the departure 
from scale invariance up to I = 1500 of 10% measured 
by Ai max (i.e. requiring Ai max < 0.1) we obtain for the 
various collapse schemes the corresponding allowed range 
of values for the parameters A and B. The results from 
these analysis are presented in the tables HI ITT1 and ITTT1 We 
see that the restriction of range in B becomes weaker for 
larger values of A, something that can be described by 
stating that the earlier the collapse occurs the larger the 
possible departures from the behavior n^k = constant. 
We note that we can recover the range of times of collapse 
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FIG. 2: Semilog plot of |ai m | 2 (Ci(fc)) for different values of (A,B), representing how robust is the scheme of collapse when it 
departs from Zk. constant. The abscissa is I until I = 2600 



for the different values of A and B. We can solve N (x) = 
A/x + B, therefore = A/k + R D B. Note that 

Rd is the comoving radii of the last scattering surface. 
Considering the radial null geodesies we find Rd = Vo — 
rjd, where rjd is the time of the decoupling. The decoupling 
of photons occurs in the matter domination epoch, so 
we can use the expression for Rd in terms of the scale 
factor, using the corresponding solution to the Friedman 



equation 

Rd = 4~ (! - V^) > ( 33 ) 

where we have normalized the scale factor so today is 
uq = 1, so, ad = a(rjd) — 10 -3 and Hq is the Hubble 
variable today. The numerical value is Rd = 5807. 31/i -1 
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FIG. 3: Semilog plot of \ai m \ 2 (C2(k)) for different values of (A,B), representing how robust is the scheme of collapse when it 
departs from Zk. constant. The abscissa is I until I = 2600 



Mpc. Henceforth 

A 2B 

l^(fc)| = ¥ + — (34) 

Thus, we can use this formula and calculate the collapse 
time of the interesting values of k we observe in the cos- 
mic microwave background (CMB), namely the range be- 
tween 10 -3 Mpc -1 < k < 1 Mpc" 1 . These modes covers 



the range of the multipoles I of interest: 1 < I < 2600, 
were we made use of the relation [24l| I = kRjy. The col- 
lapse times for these modes can be regarded as the times 
in which inhomogencities and anisotropics first emerged 
at the corresponding scales. These collapse times are 
shown in figure ([8]) for the best values of (A, B) given 
in the tablesfH ITT1 HID) . 

We can compare the value of the scale factor at the 
collapse time a(ry^), with the traditional scale factor at 
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FIG. 4: Semilog plot of \ai m \ 2 (C W i gner (k)) for different values of (A,B), representing how robust is the scheme of collapse when 
it departs from z k constant. The abscissa is I until I = 2600. 



"horizon crossing" that marks the "quantum to classi- 
cal transition" in the standard explanation of inflation: 
af: . The "horizon crossing" occurs when the length cor- 
responding to the mode k has the same size that the 
"Hubble Radius", HJ 1 , (in comoving modes k = aHj) 
therefore, = a(r]^) = = . Thus the ratio of 
the value of scale factor at horizon crossing for mode k 



and its value at collapse time for the same mode is 
n H 

-^ = krj c k (k) = A + BR D k = A + Bl. (35) 

ab- 
using the best-fit values for the different collapse 
schemes, we can plot the e-folds elapsed between the 
modes collapse and its horizon crossing. As we can see in 
the figure ^ this quantity changes -at most- of one order 
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(d)Ci(fc), A = 10 




(e)Ci(fc), A = 1000 

FIG. 5: Plot showing how the integral of |ai m | 2 (Ci) varies respect changes in B (10~ 4 — 10), keeping A fixed. Both axes B and 
I are in logscale. See the main text for a more extensive explanation. 



of magnitude in the range for k for the values of A and B 
that were considered more reasonable, i.e. > a£, the 
time of collapse r)i ~ 10~ 3 rj^ in this range. The door is 
clearly open for a more in detailed analysis and compar- 
ison to the the actual empirical data, whereby one could 
hope to extract robust information of the type discussed 
above. 



V. DISCUSSION 

We have considered various, relatively ad hoc recipes for 
the form of the state of the quantum inflationary field, 
that results, presumably from a gravitationally induced, 
collapse of the wave function. The breakdown of unitarity 
that this entails, is thought to be associated with dras- 
tic departures from standard quantum mechanics once 
the fundamental quantum gravity phenomena come into 
play. We have not discussed at any length this issue here 
and have focused in the present treatment as purely phe- 



12 




(c)C 2 (fc), A = 1 (d)C 2 (fc), A = 10 




(e)C 2 (fe), A = 1000 

FIG. 6: Plot showing how the integral of |a( m | 2 (C2) varies respect changes in B (10~ 4 — 10), keeping A fixed. Both axes B and 
I are in logscale. See the main text for a more extensive explanation. 



nomcnological aspects of the problem. 

The analysis of the signatures of the different schemes of 
collapse illustrate various generic points worth mention- 
ing: First, that, depending on the details of the collapse 
scheme and its parameters, there can be substantial de- 
partures in the resulting power spectrum, from the stan- 
dard scale invariant spectrum usually expected to be a 
generic prediction from inflation. Of course it is known 
that there exist other ways to generate modifications in 
the predicted spectrum, such as considering departures 
from slow roll and modifications of the inflaton poten- 



tial and so forth. In the approach we have been follow- 
ing the modifications arise from the details of a quan- 
tum collapse mechanism, a feature tied to a dramatic 
departure from the standard unitary evolution of quan- 
tum of physics that we have argued must be invoked if 
we are to have a satisfactory understanding of the emer- 
gence of structure from quantum fluctuations. In fact, 
by fitting the predicted and observational spectra, these 
sort of modifications are possible sources of clues about 
what exactly is the physics behind the quantum mechan- 
ical collapse or whatever replaces it. We saw that generi- 
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in' id" 

(&)C W ig ner (fc), A = 10 ^ 




(b)C r lfJ i ff7 i er (fc), A — 10 





(d)C*„ ig „ eT .(fc), ,4 = 10 




(e)C w ig ner (k), a = 1000 



FIG. 7: Plot showing how the integral of |a; m | 2 (Cwigner) varies respect changes in B (10 
B and I are in logscale. See the main text for a more extensive explanation. 



10), keeping A fixed. Both axes 



cally one recovers the standard scale invariant Harrison- 
Zeldovich spectrum if the collapse time (conformal time) 
of the modes is such that rj^k = constant J26[. On the 
other hand and as shown in detail in Q the simple gen- 
eralization of the ideas of Penrose about the conditions 
that would trigger the quantum gravity induced collapse 
leads precisely to the such prediction for rfe. We should 
however keep in mind that, even if something of that sort 
is operating, the stochastic nature of any sort of quantum 
mechanical collapse leads us to expect that such pattern 
would not be followed with arbitrarily high precision. In 



this regard we have studied the robustness of the various 
schemes in leading to an almost scale invariant spectrum. 
To this end we have considered in this work, the simplest 
(linear) deviations from the behavior of n^ as a function 
of k i.e. we have explored in the three existing collapse 
schemes the effects of having a time of collapse given by 
T}% = A/k+BRp. The results of these studies are summa- 
rized in figures [H El 0] and tables HI ITT1 and ITTTT so here we 
will only point out one of the most salient features: We 
note that the different collapse schemes lead to different 
types of departures of the spectrum from the scale invari- 
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(a)Ci, with A = 10, B = 1. 
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k (Mpc 1 ) 



(b)C 2 , with A = 1000, B = 1. 




k (Mpc 1 ) 



(c)C w i gner , with 7l = 0.01, -B = 10. Note how in this scheme almost all 
the modes must collapse at the same time. 

FIG. 8: Logarithmic plot in both axes of the collapse times \r]^\ (in seconds), for the three schemes, taking in account only the 
best values of (A,B) in the range of 10~ 3 Mpc -1 < k < 1 Mpc -1 . For these plots: h = 0.7. 



ant one, for instance the schemes C2(fc) and C W i gner (k) 
lead naturally to a turning down of the spectrum as we 
increase I. 

It is worth noting that a turning down in the spectrum 
is observed in the CMB data Q, which is attributed as 
a whole, in literature to the Damping Effect p7j , i.e. to 
the fact that inhomogencities are dampened do the non 
zero mean- free-path of photons at that time of decoupling 



[151 ] . As observed in the figures ([31 0|) for some values of 
(A, B) we obtain an additional source of "damping" due 
to fluctuations in the time of collapse about the pattern 
characterized by rfj.k = constant. It is expected that the 
PLANCK probe will provide more information on the 
spectrum for large values of I, so hopefully this character- 
istic of our analysis could be analyzed and distinguished 
from the standard damping in order to obtain interesting 
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(a)Ci, with A = 10,B = 1. 
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(b)C 2 , with A = 1000, B = 1. 




10" 3 10~ 2 10~ 1 10° 
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(c)C wigne r, with A = 0.01, B = 10. 



FIG. 9: Semi Logarithmic plot of the number of e-foldings between a k and a£ for the three schemes, taking in account only 
the best values of (A, B) in the range of 10" 3 Mpc -1 < k < 1 Mpc -1 . For these plots h = 0.7. 



constraints on the parameters (A, B). In fact we believe 
that one should be able to disentangle the two effects, be- 
cause in the cases in which our model leads to additional 
damping in the spectrum, it also predicts that should be 
a rebound at even higher values of I (see figures [2J [HI HI) - 

However, the most remarkable conclusion, illustrated by 
the present analysis, is that by focussing on issues that 
could be thought to be only philosophical and of prin- 
ciple, we have been lead to the possibility of addressing 



issues pertaining to some novel aspects of physics which 
could be confronted with empirical observations. Fur- 
ther and more detailed analysis based on direct compar- 
isons with observations are indeed possible, and should be 
carried out. This together with the foreseeable improve- 
ments in the empirical data on the spectrum, particularly 
in the large I region, and the large scale matter distribu- 
tion studies, should permit even more detailed analysis 
of the novel aspects of physics that we believe are behind 
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TABLE I: Robustness of Ci when the parameters (A, B) were 
varied from 1CT 4 < A < 10 3 and 1(T 4 < B < 10 



the origin of structure in our universe. 
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Further issues about the identification of ensemble aver- 
ages with time averages, which are valid under the er- 
godic hypothesis can be raised in the cosmological con- 
text where assumptions of ergodicity and equilibrium 
seem to be much less justified. 

Inflation could work if we don't start -strictly- with this 
condition, but after some e-foldings the universe reaches 
this stage. 

We find this wording unfortunate because it leads peo- 
ple to think that something is fluctuating in the sense of 
Brownian motion, while a wording such as "quantum un- 
certainties" would evoke something like the wave packet 
associated with a ground state of an harmonic oscillator 
which is a closer analogy with what we have at hand. 
As we will show, the relevant quantities that one is inter- 
ested in computing are determined once one characterizes 
the time of collapse and describes the state after the col- 
lapse in terms of the expectation values of the field and 
momentum conjugate variables. 

At this point, in fact, all we require is the specification 
of the expectation values of certain operators in this new 
quantum state. 

This might seem awkward to some readers. It is worth 
then emphasizing that our view is that, in contrast with 
what happens with other fields, the fundamental degrees 
of freedom of gravitation are not related to the metric 
degrees of freedom in any simple way, but instead the 
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latter appear as effective degrees of freedom of a non- 
quantum effective theory. Therefore, the quantum un- 
certainties (we feel "uncertainties" is a more appropri- 
ate word than "fluctuations", as the latter suggest that 
something is actually changing constantly in a random 
way) associated with the gravitational degrees of free- 
dom are most naturally thought as not having a metric 
description (as occurs for instance in the Loop Quan- 
tum Gravity program where the fundamental degrees of 
freedom are holonomies and fluxes), and thus that the 
metric can appear only at the classical level of descrip- 
tion, where it satisfies something close to the semiclassi- 
cal Einstein equation. In other words, from our point of 
view, it would be incorrect to think of the quantum un- 
certainties of the metric as appropriate description of the 
quantum aspects of gravitation, and much less, as satis- 
fying Einstein's equations. From our point of view, this 
would be analogous to imagining the quantum indetermi- 
nacies associated with the ground state of the hydrogen 
atom, as described in terms a perturbation of the orbit of 
an electron in Hydrogen atom, and satisfying Keppler's 
equations for the classical Coulomb potential. For more 
details about these point of view see [ES^I- The reader 
should be aware that this is not a view shared by most 
cosmologists. 

This is just a mathematical evaluation device and no as- 
sumption regarding the existence of such ensemble of uni- 
verses is made or needed. These aspects of our discussion 
can be regarded as related to the so called cosmic vari- 
ance problem. 

Note here again the difference with the standard treat- 
ment of this part of the calculation which calls for the 



evaluation of the expectation value {8<f>kS<f>k')* on the 
vacuum state which as already emphasized is completely 
homogeneous and isotropic. 

[24] The relation between the angular scale and the mul- 
tipole I is 9 ~ n/l. The comoving angular distance, 
dA, from us to an object of physical linear size L, is 
d,A = LI (off). L/a ~ 1/fc, (1a = Rd if the object is in 
the LSS, and using the first expression in this footnote, 
we get I — kRo- 

[25] The reader should keep in mind that our parametrization 
of the inflationary regime has the conformal time running 
from large negative values to small negative values 

[26] This resembles the condition that is sometimes consid- 
ered in the context of the so called trans-plankian prob- 
lem. There is however an important difference of what is 
supposed to occur at the (conformal) time that appears 
in this condition. In addressing, the trans-plankian prob- 
lem the time indicates when the mode actually comes 
into existence. In contrast , in our approach, the mode 
has existed always -modes are not created or destroyed- 
, but the state of the field in the corresponding mode 
changes (or jumps) from the Adiabatic vacuum before 
the condition to the so called post-collapsed state after 
this, or a similar condition, is reached. 

[27] This effect basically is a damping for the photon den- 
sity and velocity at scale k at the time of decoupling 
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by a factor of e D , where ko is the diffusion scale 
and depends in the physics of the collisions between elec- 
trons and photons. Accordingly, the C; spectrum is also 
damped as e~ l ^ D where Id ~ kodA^d) ~ 1500, for 
typical cosmological parameters. 



